Low-energy Effective Theory for One-dimensional Lattice Bosons near Integer Filling 
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A low-energy effective theory for interacting bosons on a one- dimensional lattice at and near 
integer fillings is proposed. It is found that two sets of bosonic phase fields are necessary in order to 
explain the complete phase diagram. Using the present effective theory, the nature of the quantum 
phase transitions among various phases can be identified. Moreover, the general condition for the 
appearance of the recently proposed Pfaffian-like state can be realized from our effective action. 
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I. INTRODUCTION 



Much effort has been devoted to understanding the ef- 
fects of competing interactions on quasi-one-dimensional 
systems. Recent advances in loading ultracold bosonic 
atoms into an optical lattice lead to the realization 
of one-dimensional (ID) lattice boson systems, and in- 
spire the investigation on many-body quantum phenom- 
ena therei m 1 ' 2 ' 3 ' 4 ' Besides the nearest-neighbor hopping t 
and the tunable on-site interaction U, a sizable nearest- 
neighbor interaction V is now within experimental reach 
by using the dipolar interaction among atoms^ These 
experiments raise the interest in creating and detecting 
exotic quantum phases in ID interacting lattice boson 
systems. 

For ID lattice bosons at integer filling, the phase di- 
agram obtained by the mean-field calculation includes 
three different phases. They are the Mott insulator (MI) 
for large U, a charge density wave (CDW) for large V, 
and a superfluid (SF) for large However, by using 
the Density Matrix Renormalization Group method, it 
is found recently that, at filling of n = 1, there exists a 
novel insulating phase lying between the CDW and MI 
phases in the weak coupling region^ The phase transition 
between this new phase and the conventional insulating 
phases (i.e., the CDW and MI phases) is found to be 
of second order. Similar to the Haldane phase of quan- 
tum spin-one chains^ this new insulating phase can be 
identified by a highly nonlocal string order parameter, 
and thus called as the Haldane insulator (HI) phase i- 
Another investigation employing quantum Monte Carlo 
simulations shows that, when t <C U < 2V, there ex- 
ists a supersolid (SS) phase lying between the n = 1 
CDW and the n = 1/2 CDW phases^ The CDW-SS 
transition is found to be of second order with dynamic 
critical exponent z = 2. Recently, due to the possibil- 
ity of a fault tolerant quantum computation based on 
non-abelian anyons, 12 people begin to search for models 
containing non-abelian anyonsji 3 - as well as for techniques 
for their detection and manipulation*^ It is pointed out 
that bosonic atoms in a ID lattice with infinite repul- 
sive three-body on-site interactions around n — 1 have a 
ground state very close to a gapless Pfaffian-like stated 



which may serve as the basis to create non-abelian anyons 
in ID systems. 



With an eye on these new developments, it is desir- 
able to have a unified understanding to the rich phases 
and the quantum phase transitions among them through 
a suitable low-energy effective theory. A common way 
of deriving low-energy effective theory for ID bosonic 
systems is to use Haldane's density-phase representation 
(also known as "phenomenological bosonization" ) 1 16 i 17 
Such an approach has been successfully used to address 
the low-energy physics of SF and the related Mott tran- 
sitions at commensurate fillings. 17 However, this effec- 
tive theory contains only one set of collective variables 
which describe long-wavelength phase and density fluc- 
tuations, and can show only one phase transition be- 
tween two phases. That is, it fails to show two phase 
transitions found in Ref. U Moreover, in this simplest 
one-component hydrodynamic effective action, the CDW 
and the SF order parameters compete (or, are dual to) 
each other. Therefore, from such an effective theory, it is 
hard to understand the origin of the SS phase reported 
in Ref. [Ill, where the CDW and the SF orders coexist. 
In the present work, a low-energy effective theory with 
two sets of bosonic phase fields is proposed. Within the 
present framework, not only the nature of the phase tran- 
sitions and the corresponding critical modes can be ex- 
plicitly identified, the phase diagrams reported in Refs. 
and nil can also be easily understood. Moreover, we find 
that the SS phase in ID has at most an algebraic long- 
range CDW order and has an algebraic boson-paired su- 
perfluid order instead of the ordinary one. Furthermore, 
the general condition for the appearance of the gapless 
Pfaffian-like state can be realized within the present ap- 
proach. This observation will be of help for searching 
gapless Pfaffian-like state in more realistic systems. We 
believe that the effective theory shown below provides a 
general description of the low-energy physics for ID in- 
teracting lattice bosons near integer filling, and it serves 
as a starting point for further investigations on the static 
and dynamical properties of these systems. 
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II. THE EXTENDED BOSON HUBBARD 
MODEL AND ITS EFFECTIVE THEORY 



tive Hamiltoniaii 22 ' 23,24 



The lattice bosons with a nearest-neighbor hopping t, 
an on-site interaction U, and a nearest-neighbor interac- 
tion V can be described by the extended Bose-Hubbard 
model (EBHM) 
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where b\ is the bosonic creation operator at site i, rii = 
b\hi is the number operator, and Srii = — n measures 
deviations of the particle number from a mean filling n. 
Since our main concern is in discussing the phase tran- 
sition among various insulating phases, where the local 
particle number fluctuations should not be strong, near 
an integer filling no (i.e., n ~ no), one can truncate the 
bosonic Fock space to only three local states with par- 
ticle numbers no and no ± l,£ i 15 i 18 i 19 i 20 Now the lattice 
boson model in Eq. (TT]) can be mapped onto a spin-one 
chain in the truncated Hilbert space, 
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1 spin operators, and the parameter 
1 is a measure of the 'particle-hole 



asymmetry'. For simplicity, £ is set to zero, which is 
valid at least when no 3> 1 . From the viewpoint of the 
bosonization theory, the main physics is indeed not mod- 
ified by the effect due to the ^-dependent termsi^l By 
further using the composite spin representation for the 
5 = 1 spin operators , 22 ' 23 where they are rewritten as a 
sum of two commuting spin-1/2 species, the problem of 
spin-one chain becomes that with two coupled spin-1/2 
chains. After applying the Jordan- Wigner transforma- 
tion and bosonization procedure for these spin-1/2 oper- 
ators, one finally arrives at the following bosonized effec- 



where $ s , 6 S and <& a , <d a are the symmetric (called as 
the charge mode) and the antisymmetric (called as the 
neutral mode) combinations of the bosonized fields for 
the two spin-1/2 operators, respectively. The bosonized 
forms of the original lattice boson operators are 




(5n,- 



\ -d x $ s - (-l) l/o - smV2n$ s cosV^V 
V 7r ira 



(5) 



with x — ja, where a is the lattice constant and plays 
the role of a short-distance cutoff. In the derivation of 
Eq. ([Jl, some irrelevant (or less relevant) terms have been 
dropped. As a result, the charge and the neutral modes 
become decoupled. We note that, while the method of 
phenomenological bosonization can give identical form of 
the charge part in Eq. ([4]),— the neutral mode is missed 
in that approach. That is the reason why the conven- 
tional theory fails to account for the complete phase di- 
agram of the interacting lattice bosons near the integer 
fillings. As shown in Refs. I22I and [23|, assuming the cou- 
plings in the two coupled spin-1/2 chains being weak, 
the values of the parameters in Eq. (T5|) can be obtained 
as follows: g\ = —52 = (2V — U)/2n 2 a, g 3 = —not/na; 
the Luttinger parameters K s = [1 + (U + %V) /imot]^ 1 / 2 , 
K a = [1 — {U — 2V)/nn$t\~ 1 / 2 ; and the sound mode ve- 
locities u s = uq/K s , u a — uo/K a with uq = note. Since 
the present effective Hamiltonian is reached under sev- 
eral approximations, the values of the above parameters 
should not be treated very seriously. Here we would like 
to consider Eq. ^ as a phenomenological theory. It pro- 
vides a universal low-energy physics of the ID interacting 
lattice boson models near integer fillings, such that the 
nature of the phase transitions and the long-distance be- 
haviors of the major correlation functions at each phase 
can be predicted. As shown below, this effective Hamil- 
tonian does explain the results reported recently ^ 11 ' 15 It 
reinforces our belief in the validity of the present frame- 
work. 
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FIG. 1: The schematic phase diagram at integer filling n = 
no in the weak-coupling regime. The phase boundaries of 
Kosterlitz-Thouless, Gaussian, and Ising types are denoted 
by the solid, dashed, and dotted lines, respectively. 





SF 


MI 


HI 


CDW 




> 1 


< 1 


< 1 


< 1 


K a 


>§ 


>h 


>! 


<h 


,9i 


irrelevant 


< 


> 


> 


92 


irrelevant 


irrelevant 


irrelevant 


< 


<*.) 


not pinned 





vT 


V 8 


<*»> 


not pinned 


not pinned 


not pinned 





<e a > 











not pinned 



TABLE I: The Luttinger parameters and the roles of various 
couplings in different phases at integer fillings. The values of 
the fields & s ,a and O a indicated in this table are the values 
at which these fields were pinned in each phase. Notice that 
throughout our discussions, we take the coupling 33 < (fol- 
lowing Ref. [22|), and it is irrelevant only in the CDW phase 
above. 



III. THE PHASES AND PHASE TRANSITIONS 
OF THE EXTENDED BOSON HUBBARD 
MODEL 

Through the discussion on the relevance of each non- 
linear term in Eq. (UJ), the phase diagram can be obtained 
as shown in Fig. [TJ which looks similar to that reported 
in Ref. in the weak-coupling regime. According to the 
work in Ref. [22|, better agreement can be reached if one 
keeps more local states (say, with particle numbers no, 
no±l, and no ±2) and maps the interacting lattice boson 
model to a spin chain with higher integer spin^S, How- 
ever, we must emphasize that one of the most important 
feature of the bosonization analysis in Ref. [13 is that 
the properties of the corresponding spin model do not 
dramatically depend on the value of the spin quantum 
number S, but only on whether S is integer or half-odd 
integer. This suggests that the main results obtained 
below from our analysis of the Hamiltonian in Eq. (|U) 
should remain unchanged even when the local particle 
number fluctuations become stronger. The phase bound- 
ary separating the SF and the MI phases is determined 
by K s = 1 . This comes from the relevance of the g\ term 
in the charge sector of Eq. (0]) , and thus gives a type of 
Kosterlitz-Thouless transition^ The MI-HI phase tran- 
sition is given by g\ = 0. It is a Gaussian transition, 
which separates two gapped phases. The HI-CDW tran- 
sition occurs when K a — 1/2. It is due to the com- 
petition of the two relevant terms (i.e., the 32 and the 
53 terms) in the neutral sector of Eq. and thus is 
a Ising-type transition^ By using the bosonized forms 
of the original lattice boson operators in Eq. |(SJ), the 
correlation functions at long distance can be calculated. 
In the SF phase, because the charge mode is gapless 
while the neutral field Q is pinned, the density fluctu- 
ation is gapless and the single-particle correlation func- 
tion shows an algebraic long range order (LRO) with an 
enhanced superfluid correlation: (b]bj) ~ \i — j| -1 / 4Ks . 
True LRO can develop when K s ^00. In the CDW 
phase, both of <& s and <I> Q are pinned. In this case, while 



8^/n t 




SS ! SF 

FIG. 2: The schematic phase diagram near integer filling n ~ 
no . The solid lines denote the commensurate-incommensurate 
transition with z = 2. 



the single-particle correlation functions decay exponen- 
tially, the staggered part of the density-density corre- 
lation shows LRO: (SriiSrij) ~ (— This justifies 
the identification of the CDW phase. In the MI and 
the HI phases, both $ s and 9 a are pinned, and both 
the single-particle and density correlation functions de- 
cay exponentially at long distance. Both in these two 
gapped phases, the lattice translation symmetry does not 
break. The only difference is that the <& s field is pinned 
at and in the MI and the HI phases, respec- 

tively. As pointed out in Ref. 0, one can distinct the HI 
phase from the MI one via a nonlocal string order pa- 
rameter defined by O ~ (5rii exp (iTtYjk=i SnkjSrij^. 
Within the present bosonization approach, such 
a string order parameter can be expressed as 
O — (sin (y/^$ s (xi)) sin (^/|$ s (xj)))rS which is com- 
pletely independent of the neutral mode. Thus it is in- 
deed nonzero in the HI phase with $ s being pinned at 
($ s ) = consistent with the findings in Ref. 0^ 

To make the relations between the relevant couplings and 
the pinned fields in each phase more explicit, we summa- 
rize the above discussions in table [I] 

Up to now, our discussion has focused on systems ex- 
actly at integer fillings. The deviation from integer fill- 
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ing can be studied by turning on a nonzero chemical po- 
tential Sfi and adding a term Sfi ^ i Srii in the Hamil- 
tonian in Eq. ([1]). (5/i is measured from the middle of 
the lobe.) In the bosonized representation, this term 

takes the form 5fi^J~^d x $ s , which affects the charge mode 

onlyi^i The d x $> s term tends to create a finite number 
of kinks in the $ s configuration, and it will compete 
with the Sine-Gordon term cos \/8it$> s which tends to 
pin down $ s . When S/i is comparable to the magnitude 
of the charge gap of the insulating phases, such a term 
will trigger a well-known commensurate-incommensurate 
transition^ As long as K s is not too small such that 
the breather mass gap is larger than the Sine-Gordon 
kink mass gap, the universal features of this transition is 
well captured by the Luther-Emergy limit of the charge 
mode, where the bosonized effective Hamiltonian can be 
further mapped to a free massive Dirac fermion theory 
with a chemical potential term^ At the tranistion point 
5 /j, = S n c such that the charge gap closes, the gapless 
charge mode excitation density obeys the scaling rela- 
tion p ~ (Sfi — (5/ic) 1 / 2 , and the dynamic critical expo- 
nent of this transition becomes z = 2. After the transi- 
tion, the charge mode is described by a Gaussian model 
with nonzero stiffness, which accounts for the SF behav- 
ior in ID. Therefore, starting from either the MI or the 
HI phase, nonzero 5(i can induce a transition to the SF 
phase. On the other hand, after the transition from the 
CDW phase, even though the gapless charge mode shows 
the SF behavior, the CDW-like behavior can remain, and 
thus the SS phase is obtained. It is because the neu- 
tral sector is unaffected and the $ a field is still pinned 
in this case even after the transition. According to this 
analysis, it is tempted to interpret the SS phase near 
integer fillings as the descendant of the CDW phase at 
integer fillings. Thus our effective theory does explain 
the occurrence of a z = 2 second-order CDW-SS phase 
transition found in Ref. [If, It is noted that, since the 
the charge mode is gapless in the SS phase, a power- 
law behavior in the density-density correlation function 
(— (SriiSrij) oc \i — j\~ Ks is expected. That is, the 
SS phase has at most an algebraic long-range CDW or- 
der. Because K s <C 1 in the SS phase, such a power-law 
decay may not be examined easily by finite-size numer- 
ical calculations. Another interesting feature of the SS 
phase is that, while the superfluid stiffness is nonzero, 
the single-particle correlation (b\bj) decays exponentially 
due to pinning of the <& a field (or the a field being dis- 
order). However, the boson-pair correlation {{b\) 2 b 2 ) in- 
deed shows an algebraic long range order. Thus the sys- 
tem behaves somewhat like a "boson-paired superfluid'' . 
This issue deserves further examinations by other ap- 
proaches. The above discussions are summarized by the 
phase diagram shown in Fig. [21 



IV. THE PFAFFIAN-LIKE STATES 



Besides the application to the ID EBHM, a recently 
proposed gapless Pfaffian-like stated as the ground state 
of three-body-hard-core bosons in ID lattice can also be 
understood within the present framework. Due to the in- 
finite three-body on-site repulsion, only three local states 
with occupation numbers n = 0, 1, 2 are necessary. Sim- 
ilar to the previous discussion, the interacting lattice bo- 
son model in Ref. [TH can again be mapped to a spin-one 
chain in Eq. © with U = V = and n = 1. That 
is, the lattice boson model studied in Ref. [H is just a 
special case of the lattice spin-one Hamiltonian studied 
in the present paper. Therefore, the proposed gapless 
Pfaffian-like state should correspond to one of the gap- 
less ground states described by the present continuum 
effective field theory. When U = V — and no = f, 
ignoring the ^-dependent terms in Eq. ^ as we did be- 
fore, the low-energy effective Hamiltonian becomes that 
in Eq. 0$ with the Luttinger parameters K s = 1 and 
K a > 1/2. Thus it indicates that the gapless Pfaffian- 
like state should be identified as the ground state laying 
on the SF-MI phase boundary, where the charge sector is 
gapless and the a field is pinned. As a nontrivial cross- 
check of this observation, we show that the correlation 
functions calculated by the Phaffian-like wave function 
under a completely different approach can indeed be re- 
produced within the present framework. By using the 
bosonized forms Eq. ([5]) of the original lattice boson op- 
erators, for the ground state laying on the SF-MI phase 
boundary, the single-particle and boson-pair correlation 
functions behave at long distance 



|i-i|- 1/4 *', ((bt) 2 b 2 j)~\i-j\- K * (6) 



with K s = 1. Notice that, on the SF-MI phase bound- 
ary, the exponent of the single-particle correlation func- 
tion is always four times larger than that of the boson- 
pair correlation function. These results agree with those 
in Ref. Ila . It supports our observation that the ground 
state laying on the SF-MI phase boundary is nothing but 
the Pfaffian-like state. Since a direct strong three-body 
interaction among atoms are rare in nature^ the au- 
thors in Ref. [H suggested to use a Feschbach-resonant 
system away from the resonance to generate an effective 
three-body interaction among atoms, thus realizing this 
novel state. According to the present discussion, we find 
that strong three-body on-site repulsion is not essential 
for the appearance of the Pfaffian-like state. Instead, 
ID Interacting lattice boson models can have a gapless 
Pfaffian-like ground state, as long as they are tuned to- 
ward the critical points of the SF-MI phase transition. 
Thus it is plausible to realize such ID Pfaffian-like states 
in an ordinary boson Hubbard model without borrowing 
the mechanism of Feshbach resonance. 
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V. CONCLUSIONS 

To summarize, a low-energy description, which suc- 
cessfully capture the rich phases of the ID lattice bo- 
son model near integer fillings, is constructed. We be- 
lieve that the present effective theory provides general 
applications to related problems. For example, paral- 
lel approaches can be used to study models of coupled 
ID hard-core lattice bosons, where similar effective theo- 
ries can be constructedi 33 ' 34 i 35 ' 36 i 37 ' 38 When local density 
fluctuations become stronger, more local states with dif- 
ferent particle numbers should be kept. Following the 



above reasoning, an effective theory with more bosonic 
phase fields should be reached. It is interesting to see if 
more novel phases can be found in this case. 
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